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Abstract 

Using 281 fb~^ of data from the Belle experiment recorded at or near the T{AS) resonance, 
we have measured the rates of the "wrong-sign" decays — > K~^tt~tt^ and K~^tt~tt~^tt~ 
relative to those of the Cabibbo-favored decays — s- K~tt'^tt^ and K~TT~^7r^ir~ . These 

wrong-sign decays proceed via a doubly Cabibbo-suppressed amplitude or via mixing; the 

latter has not yet been observed. We obtain i?ws(i^vr7r°) = [0.229 ± 0.015 (stat.)j:[3:|]J^ (sys.)]% 
and RwsiK^-^) = [0.320 ± 0.018 (stat.)l:[|:[]J| (sys.)]%. The CP asymmetries are measured to be 
-0.006 ± 0.053 and -0.018 ± 0.044 for the K+tt'tt^ and J^+vr-vr+vr- final states, respectively. 

PACS numbers: 12.15.Ff,13.25.Ft,14.40.Lb 
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Studies of mixing in the K^-K^ and B^-B'^ meson systems |l| have had an important 
impact on the development of the Standard Model (SM). The latter allowed the top quark 
mass to be predicted prior to its direct observation. In contrast, the D^-D^ mixing rate 
is strongly suppressed by Cabibbo-Kobayashi-Maskawa (CKM) factors and the GIM mech- 
anism [2]; the SM predicted rate is far below current experimental upper limits. Obser- 
vation of mixing significantly larger than this prediction could indicate new physics 
Previously, D^-D'^ mixing has been searched for in "wrong-sign" (WS) — >• i^+vr" de- 
cays 0, B 3, in WS D° ^ K+TT-TT^ and D° ^ K+ti-ti+ti- decays 0, 0, and in Dalitz- 
plot analyses of — >^ K^vr+vr" decays 0]. Here we investigate the WS multi-body modes 

^ K'^TT" (nil) with a data sample more than 30 times larger than that of previous 
studies. These modes can arise from a mixing into and subsequently decaying via 
the "right-sign" (RS) Cabibbo-favored (CF) decay ^ K^7r~{n7r). The final states can 
also arise from a doubly-Cabibbo-suppressed (DCS) amplitude; the ratio of DCS decays to 



CF decays can be used to measure the CKM phase 03 in B^ D^K'^ [It 

In this Letter we present measurements of the ratio of rates for WS to RS de- 
cays, i?i^s^"""°^ = r(DO ^ K+7r-7rO)/r(DO i^-TT+vrO) and i^^^g^"""^"") = r{D^ ^ 
K~^7T~7r~^'!T~)/r{D^ ^ K^n^-n-^Ti^). Assuming negligible CP violation, this ratio is given 
by 

i?ws = Rd + 4^Dy' + \{^'^ + y'^) , (1) 

where is the ratio of the magnitudes squared of the DCS to CF amplitudes; and x' 
and y' are "rotated" versions of the mixing parameters x = Am/F and y = AF/2F: x' = 
X cos S + y sin 6 and y' = y cos 6 — x sin 6, where 6 is an effective strong phase difference 
between the DCS and CF amplitudes The parameters x and y are mode-independent, 
depending only on the differences in mass (Am) and decay width (AF) between the two 
D^-D^ mass eigenstates, and their mean decay width (F). 

The data sample consists of 281 fb~^ recorded by the Belle experiment at KEKB an 



asymmetric e+e~ collider operating at or near the T(4S') resonance. The Belle detector is a 
large-solid-angle magnetic spectrometer consisting of a silicon vertex detector (SVD), a 50- 
layer central drift chamber (CDC), an array of aerogel threshold Cerenkov counters (ACC), a 
barrel- like arrangement of time-of-flight scintillation counters (TOF), and an electromagnetic 
calorimeter (ECL), all located inside a superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return outside the coil is instrumented to detect Kr mesons 
and to identify muons (KLM). The detector is described in detail elsewhere jislllq. 

We consider the decay chain D*~^ D'^nf K7r{nn)7rf , where the "slow" pion 7r+ has 
a characteristic soft momentum spectrum. The charge of vr^ is used to identify whether a 
OT was initially produced. We require that all tracks have at least two SVD hits 
in both r-0 and z coordinates. We use information from the TOF, ACC, and CDC to 
select kaons (pions) with momentum dependent efficiencies of 80-95% (90-95%) and pion 
(kaon) misidentification probabilities of 5-20% (15-20%). To suppress background from 
semileptonic decays, we remove tracks identified as electrons (muons) based on ECL (KLM) 
information. We select 7r° candidates that satisfy 118 MeV/c^ < M^^ < 150 MeV/c^ (±3cr 
in resolution); we then apply a mass constrained fit for the photons. We require photon 
energies to be larger than 60 (120) MeV in the barrel (endcap) region. 

^ K~^7r~7T^ candidates are reconstructed by combining two oppositely-charged tracks 
with a 7r° candidate having p > 310 MeV/c in the center-of-mass (CM) frame. The K'^n'ir'^ 
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invariant mass is required to be in the range 1.78-1.92 GeV/c^ (±60" in resolution). To 
reject background from — > K'tt^h^ in which the K is misidentified as vr and the vr as 
we calculate 771x7^0 with the K and vr assignments swapped and reject events having 
"^i^TTTrO (swapped) m the range 1.78-1.90 GeV/c^ 

^ K^7T^7i^7T~ candidates are formed from combinations of four charged tracks; ttLj^^^ 
is required to be in the range 1.81-1.91 GeV/c^ (±7o"). To reject background due to misiden- 
tification of D° — > K~n~^7r~7r~^ , we calculate ttlj^^^ with the K and vr assignments swapped 
and reject events satisfying |?^i^37r(swappcd) ~ ''^doI < 20 MeV/c^. The Cabibbo-suppressed 
decay D'^ ^ K^K~^7i~ followed by K^^vr'^'Tr" can also mimic the WS signal; to reject this 
background, we calculate m^+^- for both oppositely-charged pion combinations and reject 
events satisfying |m^+7r- — m^^(,| < 16 MeV/c^. 

The charged daughters are required to originate from a common vertex. The 
momentum vector is extrapolated back to the interaction point (IP) profile and a production 
vertex is determined. The D*~^ candidate is then formed by combining the candidate 
with a vr^. We refit the vr^ track, requiring that it intersect the production point; 
this greatly suppresses combinatorial background and improves the resolution on the energy 
released in the D* decay, Q = ^^^^--^ - M k+-k- {mz) - fn^+ . For D* + ^D°7r+ decays, Q 

is only 5.85 MeV (slightly above threshold) and provides substantial background rejection. 
We subsequently require Q<12 MeV, which is > 99% efficient. 

To eliminate D mesons produced in BB events and further suppress combinatorial back- 
ground, the reconstructed D*^ momentum in the CM frame is required to be greater than 
2.5 GeV/c. Finally, we require that the P^r degree of freedom (d.o.f) resulting from the 

vertex fit, the IP vertex fit, and the vr^ track refit be satisfactory. The fraction of events 
containing multiple signal candidates is less than 3% for both modes (and is the same for 
RS and WS decays); multiple signal candidates are retained for subsequent analysis. 

We determine the RS and WS signal yields by performing binned maximum likelihood fits 
in M-Q space with M = MKTT{nn)- The signal and background distributions are determined 
using a large Monte Carlo (MC) sample [l7|. The backgrounds can be divided into three 
categories: (a) "random vr^" background, in which a random tt^ is combined with a true 
^ K^T[~ {nir) decay; (6) charm decay background other than (a); and (c) background 
from continuum e^e~ ^uu,dd, or ss production. 

The RS signal shape as predicted by MC simulation is parameterized in M with a sum 
of a double Gaussian and a double bifurcated Gaussian with common mean, and in Q with 
a bifurcated Student's t function. Background distributions are parameterized with similar 
empirical expressions determined from MC simulation. In the RS sample fit, the mean and 
width of the signal distribution are left free to vary, while other parameters are fixed to 
MC values. The relative normalizations of individual background categories are fixed to 
MC values for the D^^K+tt-ti^ fit, and left free for the ^ K+iT~7r+7r- fit. In the WS 
sample fit, the mean and width of the signal are fixed to the values obtained from the RS 
fit; the normalizations of the backgrounds are left free to vary. 

The RS sample fit obtains a signal yield of (8.683 ± 0.002) x 10^ for D'^ K~Tr+7T^ and 
(5.259 ± 0.002) X 10^ for ^ K~7T+7r+7T~ . The WS fit finds 1978 ± 104 for D^^K+Tr~7i^ 
and 1721 ± 75 for K~^it~7t^h^ . The fit results are projected onto the M and Q 

distributions in Fig. Q for K^tt^tt^ and in Fig. |21for K^7r~7r^7r~ . The hatched 

histograms show the fit results and the points with error bars show the data. 

In — > K7r{mT) decays, intermediate resonances dominate the decay rate and cause a 
nonuniform event distribution in phase space. Since RS and WS decays may have different 
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FIG. 1: Results of the M-Q fit for 1)° K+tt^tt'^, in projections onto (a) RS M^^^o with 
MeV< Q < 12.0 MeV; (b) RS Q with 1.780 GeV/c^ < M^^nO < 1.920 GeV/c^; (c) WS M^^^o 
with 5.31 MeV <Q< 6.42 MeV; and (d) WS Q with 1.844 GeV/c^ < Mj^^^o < 1.887 GeV/c^. 



resonant substructure, their acceptances may differ. We correct the event yields for accep- 
tance and reconstruction efficiency as follows. For ^ K'^tt^tt^ , we determine efficiencies 
using MC simulation in bins of [M]^^, M^^o); for ^ K'^n^'n-'^ir" , we use bins in a five- 
dimensional space comprised of the invariant mass squared for various K, ir combinations. 
We then calculate efficiency-corrected signal yields in each bin for the RS and WS sam- 
ples. The background is taken to be the overall background yield multiplied by the fraction 
falling in that bin; the distribution of background among the bins is taken from the sideband 
\Q — 5.85 MeV| > 2.0 MeV. The resulting signal yields are summed over all bins, and the 
ratio of the total signal yields gives 

i?ws- The results are R^l"" '^^ = (2.29 ± 0.15) x lO'^ 
and i^^s^"''^^" = (3.20 ± 0.18) x 10"^ where the errors are statistical only. 

The average efficiency for a mode is obtained by dividing the signal yield from the M-Q 
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FIG. 2: Results of the M-Q fit for —>■ K'^tt 7r"^7r , in projections onto (a) RS Mkz-k with 
MeV < Q < 12.0 MeV; (b) RS Q with 1.810 GeV/c^ < M^s^ < 1.910 GeV/c^; (c) WS Mkz-, 
with 5.47MeV <Q< 6.28MeV; and (d) WS Q with 1.852GeV/c2 < M^s^ < 1.878 GeV/c^. 



fit by the total efficiency-corrected signal yield; the ratio of average efficiencies (^^g) / (f^ws) 
is 1.01 ±0.05 for D^^K^-k^tx^ and 0.98 ± 0.04 for ^ K^tt^tt+tj- . 

Contributions to the systematic uncertainty on R-^s are listed in Table IH the size of 
each term is assessed by varying the analysis as described below and repeating the fits. 
Many effects cancel in the ratio due to the similar kinematics of the RS and WS modes; 
one distinction is the significant background contribution to the WS sample. We vary 
the selection criteria over reasonable ranges (the WS yield changes by ~10%); the largest 
positive and negative variations in -R-yvs assigned as systematic errors. We check the 
parameterization of the signal shape by varying the means and widths in M and Q by ila. 
We check background fractions and parameterizations by varying individual fractions and 
distribution parameters by ila; we also try alternative functional forms. We investigate 
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TABLE I: Systematic uncertainties for -R^g, in 


percentage. 


O 

Dource 


U — > K 


^TT IT 




Selection criteria 


+b.22 


-2.38 


+5.25 


-3.78 


Signal shape param. 


+0.09 


-0.10 


+0.10 


-0.10 


Background fraction 


+0.00 


-0.07 


+0.01 


-0.01 


Background param. 


+0.42 


-2.89 


+0.34 


-0.59 


Possible fit bias 


+2.23 


-0.94 


+0.91 


-0.88 


Total 


+5.7 


-3.9 


+5.4 


-4.0 



possible fit bias by fitting a large MC RS sample; the small difference between the fitted 
yield and the true number of RS events is taken as an additional systematic error. The total 
systematic error is obtained by combining the individual terms in quadrature. 

Assuming a value for x', Eq. (P) can be used to constrain i?^ as a function of y'. This 
constraint is shown in Fig. 01 for x' = and \x'\ = 0.028; the latter value is the 95% CL 
upper limit on \x'\ obtained from our previous analysis of K^n^ decays j^. Values 

of {x',y') for different decay modes would be equivalent if the strong phase differences (S) 
for the modes were equal. In the absence of mixing (i.e., x = y = 0), our measurements give 



tan^ 9, 



c 



R^iKTm"^) = (0.85 tan^ and i?^(K37r 
angle), consistent with theoretical expectations [1 

By separately fitting the and samples, we measure the CP asymmetry 



is the Cabibbo 
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WS 
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WS 
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We obtain Acp{K7nr°) = -0.006 ± 0.053 and Acp{K37t) = -0.018 ± 0.044, which are 
both consistent with zero. The systematic uncertainties are <0.01 (much smaller than the 
statistical errors) and are neglected. The first value represents a large improvement over the 
previously-published result [7|; the second value has not been previously measured. 

In summary, using 281 fb~^ of data we measure the ratio of WS to RS decay rates for 



■K'^n^TT^ and Z}°— > A'^tt^tt+tt to be 



pK + TT- 

-"-ws 



2.29 ± 0.15(stat)+J;j9 (syst) 



X 10" 



R 



K+i 
WS 



3.20 ± 0.18 (stat)+J;;3 (syst) 



X 10" 



These results are much more precise than previously-published results 0, 0, @|- The CP 
asymmetries measured are consistent with zero. 
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FIG. 3: 68.3% CL bands for i?^, as a function of y' for x' = and \x'\ = 0.028. The latter value 
is the upper limit obtained from our analysis of ^ K~^7r~ decays assuming no CP violation 0]. 
The point with la error bars is the result from the — > K~^7r~ analysis for x' = (and no CP 
violation). Note that 6 and thus x', y' may differ for the three modes. 
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